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We review some recent developments in the subject of quantum corrections to soli- 
ton mass and central charge. We consider in particular approaches which use local 
densities for these corrections, as first discussed by Hidenaga Yamagishi. We then 
consider dimensional regularization of the supersymmetric kink in 1+1 dimensions 
and an extension of this method to a 2+1-dimensional gauge theory with supersym- 
metric abelian Higgs vortices as the solitons. 
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I. INTRODUCTION: LOCAL DENSITIES FOR QUANTUM ENERGY AND 

CENTRAL CHARGE 

A characteristic theme in the work of Hidenaga Yamagishi is the exploration of quantum 
field contributions to the structure of solitons. The first effort was his masterful study of 
the effect on the Dirac electron vacuum of the choice of chiral boundary condition on the 
electron wave function at the location of a Dirac monopole [1] . This problem can be viewed 
as an appropriate limit of the case of Dirac isospinor fermions interacting with an 't Hooft- 
Polyakov monopole, in which one member of the isospinor doublet becomes extremely light. 
He showed that there is a fractional electric charge at the monopole (i.e., dyon charge), in 
the form of a polarization charge of the electron vacuum. Thus the boundary condition in 
this context expresses the vacuum angle earlier shown by Witten to imply fractional dyon 
charge in a theory without fermions [2]. 1 

An equally remarkable and original contribution was Yamagishi's introduction, in the 
context of solitons in one space and one time dimension, of the concept of local quantum 
energy density [3]. He was led to investigate this subject by the fact that different calculations 
of the quantum energy for a soliton in Af = 1 supersymmetric theory gave different results, 
both zero and nonzero. He realized that there could be energies associated with the choice of 
boundary conditions on the wave functions for field oscillations, governing their behavior at 
the edges of a large region containing the soliton. Computing local densities would eliminate 
this ambiguity, because those contingent boundary energies would be localized far from the 
soliton center. 

Yamagishi recognized the necessity of proper renormalization of the quantum perturba- 
tion theory for obtaining correct results, and described how to do this renormalization and 
make the resulting calculation systematic. His conclusion was that there is a nontrivial local 
energy density in the examples he studied, but that the integrated energy density vanishes. 
This conclusion had considerable appeal, because it implemented an identity between the 
energy and the central charge (whose quantum correction was expected to vanish [4, 5]). 
Indeed Yamagishi calculated also the local central charge density, verifying the equality of 
the two integrated densities. The aforementioned equality arises because in such theories, of 
the two supersymmetry generators (corresponding to the two components of a fermion wave 
function), one remains unbroken even in the presence of the soliton. The difference between 
energy and central charge is proportional to the norm of a state generated by action of a 
supersymmetry operator on the soliton ground state. Thus, if one such operator is unbroken, 
i.e., annihilates the ground state, the equality holds. 

Well over a decade passed after this work, during which there was little further develop- 
ment towards consensus on the correct value of the quantum mass correction for solitons in 
such 2D supersymmetric theories. The next stage began when Rebhan and van Nieuwen- 
huizen [RV] realized that a naive energy cutoff used in much of the susy kink literature was 
inconsistent with other, more carefully defined regularization methods. They used mode 
number regularization (described below) to compute the total quantum correction for a 
supersymmetric kink with periodic boundary conditions on the fermion as well as boson 



Yamagishi's analysis contained as a special case the known fact that in the presence of strictly massless 
fermions chiral rotation undoes the effect of vacuum angle, so that the latter becomes unobservablc. 
Also, in light of the main topic of the present work, it is worth noting that Yamagishi computed the 
(logarithmically divergent) energy associated with the boundary condition. 
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fluctuation wave functions [6]. 2 Of course, this nonzero result, differing from some prior 
results, was open to the objection that it might include a boundary energy contribution. 
Not long after, the same two authors joined by Nastase and Stephanov [NRSV] imposed 
'natural' or 'invisible' boundary conditions on the fermion wave functions both in the trivial 
or vacuum background and in the kink background (so that there could be no extra energy 
localized near the boundary), and reproduced one of the finite answers from the earlier pe- 
riod [8]. Graham and Jaffe [GJ] did a calculation in terms of scattering phase shifts, with 
no explicit dependence on boundary conditions, and obtained the same result [9]. 

While this new activity agreed with the result first found by Schonfeld in his consideration 
of a kink-antikink system with the two objects well separated [10], a complete understanding 
required reconciliation of the global result with the local calculation of Yamagishi. Shifman, 
Vainshtein, and Voloshin [SVV] [11], stimulated by the work of the Stony Brook-Vienna 
group, decided to tackle the question head-on by returning to the calculation of the local 
central charge density pioneered by Yamagishi. They identified an ambiguity hidden in the 
earlier work. Yamagishi had been explicit in his computation of renormalization effects, but 
had not discussed explicitly the other V required for a one loop calculation: regularization. 
The point here, which Yamagishi clearly understood, is that in calculating an energy sum 
one must subtract from the sum in the presence of the soliton the equally (quadratically) 
divergent sum in a flat, or vacuum background. 3 The subtracted quantity is only logarith- 
mically divergent, and that is the divergence compensated by renormalization. Given the 
regularization scheme, one still has to fix the finite parts in the renormalization by suit- 
able renormalization conditions. Lacking an explicit regularization scheme (possibly he had 
simple energy cutoff in mind), Yamagishi was not in a position to determine uniquely the 
finite part of the renormalization. By implementing higher- derivative regularization for the 
central charge density, SVV were able to accomplish this while explicitly maintaining super- 
symmetry. They obtained a nonzero integrated correction to the central charge, precisely 
equal to the energy shift previously calculated by Schonfeld, as well as NRSV and GJ. 

Thus the equality of mass and central charge which Yamagishi had found was preserved, 
but both quantities were shifted from zero by an extra finite piece which was not expected. 
NRSV, on the basis of their nonzero result for the quantum correction to the energy, sus- 
pected the existence of an anomaly, given the naive expectation that the quantum correction 
to the central charge should vanish. GJ did not discuss anomalies at all, but did notice that 
the terms corresponding to each state in the unregulated, unrenormalized energy-density 
sum are equal one to one with the corresponding terms in the central charge density. They 
accomplished regulation and renormalization for energy by subtracting from the divergent 
sum the divergent contribution from the first Born approximation to the phase shift, and 
in the case of fermions also the divergent part of the second Born approximation. They 
justified this choice by noting that the no-tadpole condition in the trivial sector requires a 
counterterm for the soliton mass which is precisely equal to the integral over momentum of 
the contribution from the subtracted divergent phase shifts. They took the momentum cut- 



2 They made this choice of boundary conditions following the original work for a kink in purely bosonic 
theory by Dashen, Hasslacher, and Neveu [DHN] [7]. Amusingly in the bosonic theory, the choice of 
boundary conditions makes no difference to the result. 

3 Actually in the supersymmetric system the quadratic divergence in the trivial sector cancels between 
bosons and fermions, but this just replaces the previous regulation problem by one of correctly matching 
the energy sums for the bosons and fermions in the kink sector. 
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off used to define the counterterm equal to the momentum cutoff in the sum over momenta 
of contributions from continuum states to the energy shift. Having already determined the 
counterterm (using naive energy cutoff), by this choice they implicitly defined a regulariza- 
tion scheme. To implement the equality of central charge with energy, they used the identical 
subtraction to define central charge as they had used for energy. This certainly was consis- 
tent, but it left open the question whether and if so where an anomalous contribution to the 
central charge might be identified, as is remarked at the end of the paper of SVV. In any 
case, what remained to be demonstrated was that one could deduce the regularization and 
renormalization of the central charge directly by analyzing its field-theoretic structure. SVV 
did exactly that, showing explicitly that there is an anomaly, which preserves supersymme- 
try by providing a shift in the central charge equal to that in the energy. Like the famous 
chiral anomaly this one is an automatic consequence of using a regularization scheme that 
preserves a relevant symmetry, in the present case supersymmetry. 

In an important respect the line of development begun by Yamagishi still was incomplete, 
because the local energy density corresponding to the established global energy and local 
central charge density was not computed directly. That was accomplished by Goldhaber, 
Litvintsev, and van Nieuwenhuizen [GLV2] [12] and independently by Wimmer [13]. They 
demonstrated by explicit checks the validity of a local regularization scheme called local 
mode regularization. This scheme makes computation of energy density of a soliton in 2D 
theory very clear and simple. Not surprisingly, it confirms directly for the supersymmetric 
kink the equality of energy density and central charge density proposed by Yamagishi and 
demonstrated by SVV. At the same time, local energy density is an interesting quantity in 
non-supersymmetric theories, where the definition of central charge is not clear. 

Even with the above results, there remains an open frontier in the computation of quan- 
tum energy densities, especially in the direction of higher dimensions. Clearly, the tree of 
local energy determination which Hidenaga Yamagishi planted continues to grow and bear 
fruit, a tribute to his insight and originality. 

In section II we focus on a development which was essential to the new understanding 
of local energy and central charge densities, namely, the careful regularization by the older, 
global methods of the total mass of a soliton, as discussed in GLV1 [14]. This discussion is fol- 
lowed by an explicit demonstration of the universal equivalence, in results for the anomalous 
or high-energy contribution to quantum energy of a field in the presence of some background 
potential, between higher- derivative regularization and local mode regularization. 

In section III and IV, we describe some recent work in (slightly) higher dimensions. We 
begin by embedding the supersymmetric kink as a domain wall in 2+1 dimensions and use 
this to set up a supersymmetry-preserving dimensional regularization scheme. This allows 
us to derive the anomaly in the central charge in a particularly transparent manner, which is 
found to be made possible by parity violation and a corresponding nonvanishing expectation 
value for the momentum operator in the extra dimension. In section IV we apply this method 
to the supersymmetric abelian vortex 2+1 dimension and determine its quantum corrections 
to mass and central charge. 
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II. SOLITONS IN 1+1 DIMENSIONS 



A. Fermion zero modes in global mode regularization 



One method of regularization, utilized by DHN in their investigation of the quantum 
correction to the mass of a kink or a sine-Gordon soliton in a purely bosonic theory, is called 
mode number regularization. There is a very concrete picture behind this scheme. Imagine 
introducing, in what initially is vacuum, some background potential which influences the 
motion of a fermion obeying the Schrodinger equation. Suppose further that at the start 
one had a free Fermi sea filled to some level, so that all possible states up to a maximum 
number N were occupied. Then, provided the boundary conditions on the wave functions at 
large positive and negative values of the spatial coordinate x prevent any probability leakage 
through the boundaries, the total number N will not change when the potential is introduced. 
Further, if the Fermi level is high enough compared to the magnitude of the potential, then 
there will be no level crossings near that level as the potential is gradually altered. Thus, just 
as before the potential was introduced, the first iV states will be filled, though the maximum 
energy may be changed slightly, and the low-energy spectrum may be altered dramatically, 
for example, to include bound states. To find the regulated quantum shift in energy one 
may subtract the sums of energies of the first iV states with and without the potential 
present. The net sum might still be divergent with N, in which case renormalization would 
be required as well, but in any case the effect of this regularization would be well defined. 

One now can adopt the same maneuver for a relativistic Bose field in the presence of a 
potential. Before the potential is introduced, one uses an energy regularization to say that 
all solutions of the wave equation up to a maximum energy e max are counted. There are N 
pairs of these solutions, each pair with one positive and one (equal) negative frequency. After 
introducing the potential, we again require N pairs of solutions, and those solutions are all 
the ones having energy e\ < e' max (iV). Thus we find a regulated, but not yet renormalized, 
energy, 



A subtlety arises in this calculation if the presence of a potential happens to introduce a 
zero-frequency solution, which would occur either accidentally or as a result of some partic- 
ular symmetry obeyed by the potential. For example, if the potential for Bose fluctuations 
arises from the presence of a nontrivial classical soliton solution of the equations of motion of 
the Bose field, then there is clearly a zero frequency deformation corresponding to uniform 
translation of the soliton. As mentioned above, solutions of (equal) positive and negative 
frequency come in pairs, and one pair defines one state. At zero frequency there is also 
one state, because the 'momentum' conjugate to the Bose field for the translation solution 
becomes the momentum of the soliton moving with some definite velocity. Each pair of 
solutions of the bosonic field equations (including any at zero frequency) corresponds to one 
term in (1). 

For fermions the situation is less familiar. Again, there are paired positive and negative 
frequency solutions, each corresponding to a term in (1). However, if there is only one zero 
energy solution of the fermionic fluctuation equations, this does not correspond to a term 




(1) 



i=i 



where the factor | 
mode. 



is the familiar coefficient in the zero-point energy contribution of each 
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in (1). The coefficient of this solution in the expansion of the Fermi field is an operator cq 
satisfying c\ = 1, and that Co appears in the expression for one of the two supersymmetry 
operators (the one that does not annihilate the soliton ground state). Several remarkable 
facts follow from these observations. First, such a wave function does not correspond to 
a zero mode of the soliton, because any mode requires two independent, noncommuting 
nilpotent creation and annihilation operators, and here one has instead one hermitean, 
idempotent operator (which may be represented as a Pauli a matrix). This immediately 
introduces a serious difficulty for the standard recipe of mode regularization, because the 
spectrum of nonnegative frequency solutions no longer is in one-to-one correspondence with 
the number of terms in (1). Secondly, because the one supersymmetry operator associated 
with the zero-frequency solution does not annihilate the soliton ground state, it follows that 
the ground state spontaneously breaks supersymmetry. 

There is another remarkable aspect of this phenomenon. The ground state not only is 
not annihilated by a supersymmetry operator, but becomes an eigenstate of that operator, 
if one decomposes the kink vacua \K) and Cq\K) into 775(1 +cq)\K) and 775(1 — c o)\K) [15]. 
This is remarkable because there is a discrete Z2 operator which takes the Fermi field into its 
negative, and which in other contexts could be used to prove that one cannot build a coherent 
superposition of a boson and a fermion. Now we find that an operator anticommuting with 
this Z2 operator leaves the ground state unchanged. Therefore two states with equal and 
opposite eigenvalues for the supersymmetry generator must be identified. In other words, 
the Hilbert space splits into two noncommunicating parts, which are gauge copies of each 
other [14]. This is a discrete analogue to the Higgs mechanism, in which a continuous local 
gauge symmetry is hidden because a scalar field in a nontrivial representation of that gauge 
symmetry has a nonzero expectation value. 

Something not considered in GLV1 is the generalization of this discussion about the Z2 
gauge symmetry to the case of many solitons, which would be relevant for sine-Gordon 
solitons if not for the kink (the focus of that study). 4 The Z 2 symmetry is not a local 
gauge symmetry, so that for a multisoliton system the relative signs of eigenvalues of the 
supersymmetry operator acting on different solitons might be significant. This may be 
related to a discovery of Moore and Read long before, in a 2+1-dimensional system which 
appears to describe the 5/2 state of the fractional quantum Hall effect [16]. They discussed 
'nonabelian statistics' for objects carrying electric charge =Fe/4. This kind of statistics seems 
to us interpretable in terms of the eigenvalues of the supersymmetry operator, which would 
reverse sign when one such object made a full circle around another. Specifically, for a 
system with n solitons separated from each other by a finite distance, charge conjugation 
symmetry guarantees [n/2] pairs of equal positive and negative frequency solutions, and for 
odd n one zero frequency solution. To each of the paired solutions corresponds a fermion 
creation operator, and hence there are 2t n//2 l nearly degenerate states of the n-soliton system. 
For infinite separation, there are n operators c which each have equal positive and negative 
eigenvalues. Superficially this might imply 2 n distinct states, but as we already have argued 
there is a gauge redundancy which reduces this number. The calculation for finite separation 
shows that the reduction must bring that number down to 2^ a l 2 ^. 

While these ideas are surprising and beautiful, one would like to have a clear idea what 
has happened to the principle of mode number regularization, now that the counting of 



Because a kink cannot be adjacent to another kink, pure multi-kink systems, with no anti-kinks, are not 
possible. 
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fermion modes suddenly has a lacuna. GLV1 found two different ways to deal with this 
question. First, following NRSV, one may use locally invisible (periodic or antiperiodic) 
boundary conditions in vacuum, and locally invisible (twisted periodic (TP) or twisted 
antiperiodic (TAP), where the twist refers to a chiral rotation of a wave function by it/ 2 at 
the one boundary with respect to that at the other) conditions in the presence of a kink. 
Now an assumption necessary to derive the equality of numbers of modes before and after 
the background is introduced no longer holds, because the twisting at the boundary allows 
'leakage' of mode number at the boundary. This is a concept that goes back to an article 
of Goldstone and Wilczek describing adiabatic flow of fermion number out of a magnetic 
monopole as the Yukawa coupling of the fermions to a Higgs field is chirally rotated [17]. 
In their discussion the leakage is out to infinity rather than into the boundary. The final 
result for localized charge is the same, whether the chiral rotation is applied to the mass, 
while the boundary condition is fixed, or applied to the boundary condition, while the mass 
remains constant (as in Yamagishi's work on monopole dyon charge [1]). The leakage in our 
example involves exactly half a unit of mode number, and this shift reconciles the count of 
fermion modes with the principle of mode number regularization. 

A second approach is to use fixed boundary conditions, which insures no leakage of 
mode number, but to average over a set of boundary conditions, so that in the average any 
contributions to the energy localized near the boundary cancel out. In this case, the mode 
number regularization is used for each set of fixed boundary conditions, but in some cases 
there are two zero-frequency solutions in the kink sector, one localized near the kink and 
one localized near the boundary, and in some cases there is one in the kink sector localized 
at the kink and one in the vacuum sector localized near the boundary. When the difference 
between the sectors is two solutions in the kink sector, this implies one (delocalized) zero 
mode in the kink sector. Thus the average over all boundary conditions amounts to half a 
zero mode excess in the kink sector. 

We see that the two different approaches agree: Effectively there is a half-mode at zero 
energy in the kink sector. This was discussed already by GJ (who analyzed half-modes at 
zero energy as well as nonlocalized 'half-bound' modes at the continuum threshold). What 
GLV1 added was a precise specification of the boundary conditions which make this notion 
well defined. The half-mode has the consequence for the mode sums used to compute the 
energy that at the cutoff one needs to include an excess of one half-mode in the sum for 
the vacuum energy compared to the sum for the kink energy. One may recognize this fact 
directly by matching modes so as to insure that no difference between the two sums linear 
in the energy will occur. Thus, as with many other examples of anomalies and associated 
regulation prescriptions, one may find the correct prescription and its implications either by 
focusing on the lowest energies (in this case zero modes), or on the highest energies. 

The half mode bears an interesting relation to the original discovery by Jackiw and Rebbi 
[18] that a soliton can polarize a charged fermion vacuum to localize a half unit of fermion 
number around the soliton. For neutral excitations, there is no directly observable charge 
or charge density that corresponds to a mode or to the probability density in its wave 
function. Nevertheless modes of nonzero energy contribute to vacuum zero point energy and 
energy density. Thus mode number and mode density are mathematically well defined and 
indirectly observable, making them also meaningful, even if not so tangible as the directly 
observable charge and charge density of a charged field. In other words, mode number 1/2 
has operational significance just as does charge 1/2. A further indication of this significance 
is additivity: Two configurations each possessing a zero frequency solution corresponding to 
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mode number 1/2, when brought within a finite distance of each other, accommodate one 
observable excitation with near-zero energy. Thus, even though it is not possible to define 
directly a "half-excitation" it is possible to count a half-mode, and to combine two of them 
to make not only a full mode but also a full excitation or fermion state. 

Although the invisible boundary conditions guarantee that the quantum correction is 
precisely the mass shift of the soliton, there remains even in this case a strong reason for 
averaging over each of the two types of condition, namely to enforce time-reversal symmetry. 
A chiral or twisted boundary condition produces different spectra for right-moving and left- 
moving waves. Time reversal symmetry interchanges TP and TAP boundary conditions, so 
that by averaging over the two one restores the symmetry. In the vacuum or trivial sector, 
periodic and antiperiodic conditions (being real) each separately produce spectra invariant 
under the symmetry. 

B. Local mode regularization from higher- derivative regularization 

SVV in their discussion emphasized the value of using a local regularization scheme 
which enforces supersymmetry, in particular, higher- derivative regularization. In GLV2 it 
was shown that a familiar scheme, point-splitting regularization, implies local mode regular- 
ization for energy density. Here we show directly that higher-derivative regularization also 
implies local mode regularization for energy density, thus verifying explicitly its consistency 
with supersymmetry. 

In higher derivative regularization, extra terms are added to the Lagrangian involving 
extra factors of the square of the spatial gradient of the field. This gives rise to an equation 
for small fluctuations of a Bose field, including a possible background 'potential' V(x), 5 

u 2 {l-d 2 JK 2 )ct> = -d 2 x <j) + m 2 <j) + V{x)<j) + d A x <t)/A 2 . (2) 

Our notation uses A as the regulator mass called M r by SVV. The effect of the higher- 
derivative term on the classical kink solution and the resulting potential V(x) for fluctuations 
tends to zero for large A, and therefore is ignored in our discussion, i.e., we use the V 
obtained from the classical solution of the unregulated equations. Here as in GLV2, we may 
use the JWKB approximation to estimate the wave functions at large uo. This allows one to 
determine with sufficient accuracy the phase space density of these wave functions both for 
the trivial or vacuum background with V = and the kink background with 

V = -3m 2 /2cosh 2 (ma;/2) , (3) 

where m is the mass for small fluctuations about the trivial vacuum. Then one may obtain 
a subtracted density which can be integrated to give the high-energy contribution to the 
regulated net local mode density, as we shall now see. 

Concentrate on a high-energy regime, with wave vectors \k\ > K ^> W, with u 2 = 
k 2 + m 2 /(l + k 2 /A 2 ), and take A >> K. Write the wave function as <f>(x) = e tkx e tf( - x \ and 
make a gradient expansion, f(x) = f^(x) + f^(x) + .... : 

2kfW\l + ^) = -V(x) . (4) 



5 For simplicity we focus here only on the pure Bose field case, but this time there is no difficulty in 
generalizing the procedure to include Fermi fields as well. 
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and 



2A:/ (1) '(l + ^) = -^/( )(l + 5^) = 



(5) 



Because is imaginary, it alters the modulus of 0, and because and V both vanish 
at spatial infinity, this equality relating first derivatives of these two functions implies an 
equality relating the functions themselves. 

We need now to compute the mode density difference between the case with and without 
nonzero V(x). Define normalization of the wave function <pk{x) so that <j) 2 is equal to unity 
far from the region of the potential. The regulated mode density before subtraction is then 



dp(x)/d\k\ = <f) k {xy 



1 - dl/K- 



(6) 



To lowest order in V/K 2 and V/A 2 we now compute the mode density difference 



d5p/d\k\ = 



V(x) 



1 



' k 2 

4 +i 



i 



* 2 



a' 2 ; 



(7) 



Here the first term in the large braces comes from the influence of V(x) on the absolute 
square of the wave functions, while the second term comes from the influence of V on the 
regulating denominator in Eq. (6). To first order in V, this denominator factor is 



[i - dl/K" 



, k 2 

A 2 



V(x) 



. k 2 ' 

A 2 



(8) 



Omitting terms of order 1/A, the integral of (7) over the range \k\ = K to \k\ = oo gives 
the extra contribution at high energies to the local mode density (number of modes per unit 
length) 

Pmod e(x)=V(x)/2TrK . (9) 

Note that this result is insensitive to A 3> K. 

In the original local mode regularization approach, this extra density was treated as all 
associated with momentum \k\ — K, so that the resulting 'anomalous' extra energy density 
resulting from the regularization is 



K 

£anom(a:) = — Pmodet^) = V {x) / AlX , 



(10) 



including a correction by a factor -2 to Eq.(57) of GLV2, Now, for the full calculation 
with higher- derivative regularization, we must compute flul^ d\k\(dSp/dk)(\k\/2), where 
the factor \k\/2 at the end is the zero-point energy corresponding to a given mode (neglecting 
quantities of 0(m 2 /k 2 )). This integral is easily evaluated, yielding a total high-energy 
contribution to the regulated (but not yet renormalized) energy density 



unrenormalized 



(*) = 



V(x) 
8tt 



A 2 ' 

2 + £n K2 



(11) 
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Here the anomaly term 2 in Eq. (11) bears exactly the same ratio to the term (compensated 
by mass renormalization) £nA 2 as found by SVV in their analysis of the supersymmetric 
kink, e.g., their Eq. (4.5). This ratio of course is trivial in local mode regularization, as is 
easy to see from the description above. However, in Eq. (11), the term 2 comes half from 
the first term in Eq. (7) and half from the second, while the logarithm clearly comes entirely 
from the first. 

Looking at this straightforward but still somewhat intricate computation, one may un- 
derstand why SVV focused on calculation of the anomaly in the central charge density, and 
did not attack directly the anomaly in the energy density. However, Eq. (11) shows that the 
result of local mode regularization for the anomaly is universally correct in higher-derivative 
regularization, and thus completes a circle: 

DHN used (global) mode regularization for the mass of the bosonic kink, but refrained 
from tackling the problem with fermions. Yamagishi introduced local energy and central 
charge density calculations, but did not fully address the role of regularization in these 
calculations. A variety of approaches using global methods managed to respect constraints 
of regulation and still circumvent the difficulty of contributions to quantum energy from 
boundary conditions on fermions. 

SVV made a complete calculation of the central charge density (both for the anomaly and 
for other contributions) but did not compute energy density directly, so that their method 
was not applicable in its original form for purely bosonic structures (or indeed any structure 
without supersymmetry) . GLV2 introduced an efficient method, local mode regularization, 
for computing the anomaly in the energy density of any system, but previous history of 
regularization schemes implies that one should check whether a new scheme preserves sym- 
metries which we want keep at the quantum level, in this case supersymmetry. The present 
calculation, verifying local mode regularization for arbitrary background potentials starting 
from higher-derivative regularization (which does preserve supersymmetry), brings us back 
to the beginning with a reliable local scheme for doing what DHN did globally. 

There is another interesting aspect of local mode regularization. The derivation from 
point-splitting regularization in GLV2 involves a shift in the real part of the effective wave 
number at a given value of k. However, the original method, and also the computation 
done here, involves an integration over k of the imaginary part of the wave number. This 
has a suggestive resemblance to dispersion relations connecting the real part of a scattering 
amplitude to an integral over the imaginary part, related by unitarity to a cross section. 

C. Why is rigorous, manifest supersymmetry unnecessary for correct calculation? 

We have seen that a completely rigorous and systematic way to assure consistency with 
supersymmetry in computing the quantum correction to the mass and central charge of 
a soliton is to insist on explicit supersymmetry at every step. Thus SVV adopted the 
approach of higher-derivative regularization in the ultraviolet and supersymmetric boundary 
conditions discretizing the spectrum to provide infrared regularization. They then computed 
a local central charge density in the vicinity of the kink, and its spatial integral in a large 
but finite region around the soliton, which gives the desired quantum correction. While this 
is the safest method, and also fairly straightforward because the central charge density has 
a weaker superficial divergence than the energy density, a number of other methods all yield 
the same result. 

It seems to us that this not only demonstrates that the result is correct, but also suggests 
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that there is a deeper principle behind the robustness of that result. The principle, we 
suggest, is cluster decomposition. In GLV1 it was deduced erroneously from an incomplete 
calculation that besides localized energy in the vicinity of the kink and possible localized 
energy near the boundary, there can be a delocalized energy, giving a finite contribution. 
They imposed the principle of cluster decomposition to sum over different boundary condi- 
tions, so that the delocalized energy violating cluster decomposition would cancel out. It 
was observed by Wimmer [13], and argued in detail in the light of the discrete C,P. and T 
symmetries by Goldhaber, Rebhan, van Nieuwenhuizen, and Wimmer [19] that such delocal- 
ized energy does not appear in mode number regularization for known boundary conditions. 
Cluster decomposition is satisfied automatically. Thus, a boundary condition which violates 
the supersymmetry preserved by the kink can alter the energy localized near the boundary, 
but cannot influence the energy localized near the soliton, i.e., the quantum correction to 
the kink mass. This means that if one uses boundary conditions which do not result in 
boundary energy, or if one averages over conditions so that there is no net boundary energy, 
then the resulting global calculation is reliable for the mass of the kink. This is true even if 
the boundary conditions do not preserve supersymmetry - that violation does not propagate 
into the region of the kink. 

As pointed out in GLV2, even methods whose consistency with supersymmetry has not 
been explicitly established, such as point-splitting regularization, give the same answer not 
only for the energy density but also for the central charge density. Thus these local densities 
are quite robust with respect to the choice of (local) regularization methods. This suggests 
that, despite the long history of discrepant calculations, one must make quite a large error 
to get the wrong answer for the local densities. For example, failure to specify carefully the 
regularization method can make the result indeterminate, and therefore most likely wrong. 
The history of incorrect global calculations draws attention in particular to regulation by 
energy cutoff, which was used in many of those calculations. As local mode regularization 
makes manifest, using a fixed and identical energy cutoff for kink and trivial backgrounds 
is incorrect. Litvintsev and van Nieuwenhuizen [20] proposed a method to repair energy 
cutoff, which is mathematically equivalent to replacing it with mode number regularization. 

III. DIMENSIONAL REGULARIZATION OF THE SUSY KINK 

Usually dimensional regularization breaks susy. However, the 1+1 dimensional J\f = 1 
susy kink can be embedded as a domain wall in 2+1 dimensions with the same field content 
while keeping M = 1 susy invariance. 

For the corresponding classically BPS saturated domain wall (a 1+1 dimensional object 
by itself), [21] has also found a nontrivial quantum correction to the energy density. In order 
to have BPS saturation at the quantum level in 2+1 dimensions, there has to be a matching 
correction to the momentum in the extra dimension which corresponds to the central charge 
of the 1+1 dimensional case. 

In this section we show that if one uses susy-preserving dimensional regularization by 
means of dimensional reduction from 2+1 dimensions, one indeed finds the required cor- 
rection to the extra momentum. Such a nonvanishing correction turns out to be possible 
because the 2+1 dimensional theory spontaneously breaks parity. 

By dimensionally reducing to 1+1 dimensions, the parity- violating contributions to the 
extra momentum turn out to provide an anomalous contribution to the central charge as 
obtained in Ref. [11], thereby giving a novel physical explanation of the latter [22]. This is 
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in line with the well-known fact that central charges of susy theories can be reinterpreted as 
"momenta" in higher dimensions. 

The latter statement has to be handled with care, though. The classical central charge 
stems entirely from the classical antisymmetric part of the energy momentum tensor of 
the 2 + 1 dimensional theory and thus would be missed by dimensional reduction if one 
were to start in 2 + 1 dimensions with the gravitational energy momentum tensor, which is 
always symmetric on-shell (in the absence of local Lorentz anomalies) and which contains 
the genuine momentum operator. However, it is the symmetric part of the 2 + 1 dimensional 
EM-tensor which gives the anomalous contribution to the quantum central charge. This 
anomalous contribution can be reduced to a surface term and is thus completely determined 
by the topology of the soliton background, independent of the precise field profile in the 
bulk. Therefore when we refer to the <p A kink in the following, this is just a special case of 
a more general situation. 

In the case of the susy kink, standard ('t Hooft-Veltman) dimensional regularization is 
seen to be compatible with susy invariance only at the expense of a spontaneous parity 
violation, which in turn allows non-vanishing quantum corrections to the extra momentum 
in one higher spatial dimension. On the other hand, as we shall see, the surface term that 
provides the classical central charge does not receive quantum corrections in dimensional 
regularization, by the same reason that led to null results previously in other schemes [5, 6, 8]. 
In dimensional regularization (by going up in the number of dimensions), the nontrivial 
anomalous quantum correction to the central charge operator is thus seen to be entirely the 
remnant of the spontaneous parity violation in the higher-dimensional theory in which a 
susy kink can be embedded by preserving minimal susy. 

Alternatively, we shall consider dimensional regularization by dimensional reduction from 
1 to 1-e spatial dimensions, which also preserves supersymmetry. In this case we show that an 
anomalous contribution to the central charge arises from the necessity to add an evanescent 
counterterm to the susy current [22]. This counterterm preserves susy but produces an 
anomaly in the conformal-susy current. We also construct the conformal central-charge 
current [22] whose divergence is proportional to the ordinary central-charge current and 
thus contains the central-charge anomaly as superpartner of the conformal-susy anomaly 
[11]. 

A. The model 

A real scalar field model in 1+1 dimensions with spontaneously broken Z 2 symmetry (</? — > 
—ip) has topologically nontrivial finite-energy solutions called "kinks" which interpolate 
between the two neighboring degenerate vacuum states, as for example = ±v. If the 
potential is of the form V(ip) = ^U 2 ((p) it has a minimal, J\f = (1,1), supersymmetric 
extension [23] 

C = -i [(d^) 2 + f%) 2 + jrfd^ + VMW*] (12) 

where ip is a Majorana spinor, ip = ifj T C with = — (7 M ) T C. We use a Majorana 
representation of the Dirac matrices with 7 = — ia 2 , 7 1 = a 3 , and C = a 2 in terms of the 
standard Pauli matrices a k so that ip = (^_) with real ip + (x,t) and ip~(x,t). (The reason 
for choosing 7 1 = a 3 , rather than 7 1 = a 1 , is that it diagonalizes the Dirac equation.) 
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The p 4 model is defined as the special case 

t%) = /f (y? 2 -^ 2 ), vl^^/x (is) 

where the Z 2 symmetry of the susy action also involves the fermions according to ip — > 
— (p,ip — > 7 5 -^ with 7 5 = 7°7 1 . A classical kink at rest at x = which interpolates between 
the two vacua is given by [24] 

Pk = vo tanh (jiqx/Vt^ . (14) 

At the quantum level we have to renormalize, and we shall employ the simplest possible 
scheme 6 which consists of putting all renormalization constants to unity except for a mass 
counter term chosen such that tadpole diagrams cancel completely in the trivial vacuum. 
So we set Z^ = = Z\ = 1 and /x 2 , = /i 2 + <5/i 2 , for which at the one-loop level and using 
dimensional regularization we find 

. 2 x r 2 a f dk d d k -% . f d d k 1 . . 

<J/r = XSv = X / - " — = A / - . , =i , (15) 

^ J (2ix) d+l k 2 + m 2 - le J (2n) d 2 [k 2 + m 2 } 1 / 2 

where m = U'(v) = y/2fi is the tree- level mass of elementary bosons and fermions, and 
h 2 — P - h 2 

The susy invariance of the model (12) under 

5y9 = # , 6ip = (0<p-U)e, (16) 

(with fxl replaced by /i 2 + 8fx 2 ) leads to the on-shell conserved Noether current 

i, = -(*»?+ i%0fr^ ( 17 ) 

and two conserved charges Q ± = J dx j 0± . 

The model (12) is equally supersymmetric in 2+1 dimensions, where we use 7 2 = a 1 , and 
the Noether current and charges are unchanged in form. The same renormalization scheme 
can be used, only the renormalization constant (15) has to be evaluated for d = 2 — e in 
place of d — 1 — e spatial dimensions. 

While classical kinks in 1+1 dimensions have finite energy (rest mass) M = m 3 /A, in 2+1 
dimensions they yield domain walls with a profile given by (14) and finite surface (string) 
tension M/L = m 3 /X. With a compact extra dimension one can use these configurations to 
form "domain strings" of finite total energy proportional to the length L of the string when 
wrapped around the extra dimension. 

The 2+1 dimensional case is different also with respect to the discrete symmetries of (12). 
In 2+1 dimensions, 7 5 = 7°7 1 7 2 = ±1 corresponding to the two inequivalent irreducible 
representations for 7 2 = icr 1 . Therefore, the sign of the fermion mass (Yukawa) term can no 
longer be reversed by ip — > ^ip and there is no longer the Z% symmetry tp — > —tp, ip — > ^ip. 

What the 2+1 dimensional model does break spontaneously instead is parity, which 
corresponds to changing the sign of one of the spatial coordinates. The Lagrangian is 



6 See [21] for a detailed discussion of more general renormalization schemes in this context. 



15 



invariant under x m — > —x' m for a given spatial index m = 1 or m = 2 together with ip — > — </? 
(which thus is a pseudoscalar) and ^ — > 7 m -?/>. Each of the trivial vacua breaks these 
invariances spontaneously, whereas a kink background in the x 1 -direction with ipk{— x 1 ) — 
—ip K (x 1 ) preserves x 1 = x reflection symmetry, but breaks it with respect to x 2 = y. 

This is to be contrasted with the 1+1 dimensional case, where parity (x 1 — > —a; 1 ) can be 
represented either by ip — > 7°-?/> and a true scalar ip — > </? or by ^ — > 7 1 ?/; and a pseudoscalar 
<£> — > — </?. The former leaves the trivial vacua invariant, and the latter the ground state of 
the kink sector. 



B. Susy algebra 

The susy algebra for the 1+1 and the 2+1 dimensional cases can both be covered by 
starting from 2+1 dimensions, the 1+1 dimensional case following from reduction by one 
spatial dimension. 

In 2+1 dimensions one obtains [25] 

{Q a , Q p } = 2i( 1 M )%P M , (M — 0, 1, 2) 

= 2t( 1 °H + 7 1 (P, + Z y ) + 7 2 (^ - Z x ))%, (18) 

where we separated off two surface terms Z m in defining 

P m = J d d xV m , V m = <pd m <p- i(^7°<9 m ^), (19) 

Z m = J d d xZ rn , Z m = U(<p)d m (p = d m W(<p) (20) 

with W(ip) = / dip U (<p). Note that the usual central charge density of the two-dimensional 
model, Z m , is obtained by dimensional reduction of the antisymmetric part of the three- 
dimensional energy momentum tensor. The local version of the susy algebra (18) is obtained 
by a susy variation (16) of the supercurrent (17) as follows 

T MN ~ Tr{ ri M 5 ] N ) = Tr( 7 M 7 JV 7 p ) d P ipU{ip) + symm. part 

~ e MNP dpipU{ip) + symm. part, (21) 

and the central charge density is then the momentum density T 02 in the reduced extra 
dimension. 

With our choice of Dirac matrices the supercharges and the superalgebra they generate 
are given by 



Q ± = J d 2 x[(p + dyip)^ + (d x ip ± U(ip))^}, (22) 

{Q ± ,Q ± } = 2(H±(Z X -P y )), {Q+,Q-} = 2(P X + Z y ) . (23) 

Having a kink profile in the x-direction, which satisfies the Bogomol'nyi equation d x ipK = 
—U(lpk), one finds that the charge Q + (corresponding to the terms in (16) with e~) leaves 
the classical topological (domain- wall) vacuum (ip = ipx, ip — 0) invariant. This corresponds 
to classical BPS saturation, since with P x — and P y = one has {Q + , Q + } = 2(H + Z x ) 
and, indeed, with a kink domain wall Z x jL d ~ x = W(+v) - W(-v) = —MjL d ~ x . 
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At the quantum level, hermiticity of Q ± and positivity of the Hilbert space norm imply 
a lower bound for the energy (density): 

(E|tf|E> > | (E|P y |E> | = |(E|(P y - Z X )|E)|, (24) 

where |£) denotes any state in the Hilbert space. This inequality is saturated when 

Q+|£> = 0. (25) 

Massive BPS states in 1+1 dimensions correspond to massless states in 2+1 dimensions, 
since with [H, P m ] = one has 

(P M P M ) = -\{{Q +2 Q- 2 - {Q\ Q-} 2 ) = (26) 

for BPS saturated states (25) with (P y ) = M for a kink domain wall with kink profile in 
the x-direction. An anti-kink domain wall has instead Q~\Yl) = 0. In both cases, half of the 
supersymmetry is spontaneously broken. To take into account that there is infinite momen- 
tum and energy unless the y-direction is compact with finite length L, one can formulate 
the above identities for energy and central charge per unit length or for energy and central 
charge densities. 

Omitting regularization the susy algebra in 1+1 dimensions is obtained from (18) simply 
by dropping P y as well CIS Zy SO that P x = P x . The term 7 2 Z a , remains, however, with 7 2 
being the nontrivial 7 5 of 1+1 dimensions. Identifying Z x with Z, the susy algebra simplifies 
to 

{Q ± ,Q ± } = 2(H±Z), {Q + ,Q-} = 2P X (27) 
and one obtains the quantum BPS bound 

(v>i#iv> > mzm (28) 

for any state BPS saturated states have Q + \ip) = or Q~\ip) = 0, corresponding to 
kink and anti-kink, respectively, and break half of the supersymmetry. 

C. Fluctuations 

In a kink (or kink domain wall) background one spatial direction is singled out and we 
choose this to be along x. The direction orthogonal to the kink direction (parallel to the 
domain wall) will be denoted by y. 

The quantum fields can then be expanded in the eigenfunctions, which are known ana- 
lytically for the (p A and sine-Gordon soliton [24] , times plane waves in the extra dimensions. 
For the bosonic fluctuations we have 

[-□ + {U' 2 + UU")]t] = (29) 

which is solved by 

71 = I £^ % (° M e ~ l( "~ ey) ^ + e<(w ^(*)) • (3°) 
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The kink eigenfunctions <fi k are normalized according to J dx\<f)\ 2 = 1 for the discrete states 
and to Dirac distributions for the continuum states according to J dx$* k $y = 2ir5(k — k'). 
The latter are deformed plane waves because there is no reflection in the case of the kink. 
The mode energies are 00 = ^Juj1 + £ 2 where u k is the energy in the 1+1-dimensional case. 

The canonical equal-time commutation relations [r)(x),r)(x')] = i5(x — x') are fulfilled 
with 

[a k ,e, a\, A ,} = 5 kk <5(£ - £'), (31) 

where for the continuum states 8 k ,k' becomes a Dirac delta function. 

For the fermionic modes which satisfy the Dirac equation [ ft + U']ip = 0, i.e. explicitly 

(d x + U')ij + + z(u + = 0, 
(d x - U')ip- + i(uj - £)ip + = 0, 

one finds 

^ = J -0^0, e~«™ ( O °) , bl(£) = b (-£). (32) 

The fermionic zero mode 7 of the susy kink turns into massless modes located on the domain 
wall, which have only one chirality, forming a Majorana-Weyl domain wall fermion [21, 26- 
28]. 8 

For the massive modes the Dirac equation relates the eigenfunctions appearing in the 
upper and the lower components of the spinors as follows: 

s k = —(d x + U')cf> k = —L={d x + U')(f> k , (33) 
uj k — t 

so that the function s k is the SUSY-quantum mechanical [29] partner of <p k and thus coincides 
with the eigen modes of the sine-Gordon model if <p k belongs to the <£> 4 -kink (hence the 
notation) [30]. With (33), their normalization is the same as that of the <p k . It is the 
relation (33) and the fact that it relates bosonic to fermionic modes, as well as different 
components of the fermionic modes to each other, which makes it possible to compute the 
one loop-corrections to energy and central charge without explicit knowledge of the mode 
functions. 

The canonical equal-time anti-commutation relations {tp a (x), ^(x 1 )} = 5 al3 5(x — x') are 
satisfied if (using that 4>- k (x) = <fi* k (x) and S- k (x) = s* k (x)) 

{b (£),bl(£')} = {b (£),b (-£')} = 5(£-£'), 

{hsA't'} = h,k'S(£-£'), (34) 



7 By a slight abuse of notation we shall always label this by a subscript 0, but this should not be confused 
with the threshold mode k — (which does not appear explicitly anywhere below). 

8 The mode with I = corresponds in 1+1 dimensions to the fermionic zero mode of the susy kink. If there 
are no other fermionic zero modes, it has to be counted as half a degree of freedom in mode regularization 
[14]. For dimensional regularization such subtleties do not play a role because the zero mode only gives 
scale-less integrals and these vanish. 
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and again the 8^ becomes a Dirac delta for the continuum states. The algebra (34) and 
the solution for the massless mode (32) show that the operator bo(£) creates right-moving 
massless states on the wall when £ is negative and annihilates them for positive momentum 
£. Thus only massless states with momentum in the positive y-direction can be created. 
Changing the representation of the gamma matrices by 7 2 — > — 7 2 , which is inequivalent 
to the original one, reverses the situation. Now only massless states with momenta in the 
negative y-direction exist. Thus depending on the representation of the Clifford algebra one 
chirality of the domain wall fermions is singled out. This is a reflection of the spontaneous 
violation of parity when embedding the susy kink as a domain wall in 2+1 dimensions. 

Notice that in (32) d can be only 2 or 1, for which £ has 1 or components, so for strictly 
d — 1 one has £ = 0. In order to have a susy-preserving dimensional regularization scheme 
by dimensional reduction, we shall start from d = 2 spatial dimensions, and then make d 
continuous and smaller than 2. 



D. Energy corrections 

Before turning to a direct calculation of the anomalous contributions to central charge 
and momentum, we derive the one-loop energy density of the susy kink (domain wall) in 
dimensional regularization. 

Expanding the Hamiltonian density of the model (12), 

H = \[tp + (Vy?) 2 + U\ip)\ + ^ V[7 V + U'{tp)\1>, (35) 
around the kink/domain wall, using ip = <£k + r), one obtains 

n = \m<PK? + u 2 ] - ^=U - d x (Urj) + 

+\W + (V# + \(U 2 )"V 2 } + ^V[7 V + U'}4, + 0(h 2 ), (36) 

where U without an explicit argument implies evaluation at tp — ipx and use of the renormal- 
ized n 2 . The first two terms on the r.h.s. are the classical energy density and the counterterm 
contribution. The terms quadratic in the fluctuations are the only ones contributing to the 
total energy. 9 They give [22] 

4/<0^|>-M 2 ). (37) 

When integrated, the first term, which is a pure surface term, cancels exactly the counterterm 
(see (15)), because 

Jm^auv)) = It/vr. =™J£^J~ = »* 2 . (38) 



The third term in (36) is of relevance when calculating the energy profile [11, 12]. 
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where we have used that U'(x = ±00) = ±m. 

In these expressions, the massless modes (which correspond to the zero mode of the 1+1 
dimensional kink) can be dropped in dimensional regularization as scale-less and thus van- 
ishing contributions, and the massive discrete modes cancel between bosons and fermions. 10 
Using the explicit form of <f>k(x), Sk(x) as given e.g. in [12], the rc-integration over the con- 
tinuous mode functions gives the required difference of spectral densities as 

/2m 
dx(\Mx)\ 2 - \s k (x)\ 2 ) = -6'(k) = -—— -, (39) 
k z + m A 

where 9(h) is the additional phase shift of the mode functions s k compared to 4>k- 

With the help of the susy-quantum mechanical relation (33) for the fermionic modes in 
the BPS background the integral (39) can also be computed without detailed knowledge 
of the mode functions [31]. Denoting the operator in (33) by A = d x + U' the fluctuation 
equation above (30) of the "bosonic" modes 4> k factorizes as 

tfAfa = oulfa, (40) 

where A^ = —d x + U' is the adjoint operator. Using (33) the spectral density (39) can be 
written as 

-6'(k) = [dx i\Mx)\ 2 - -^(A^YtA^)} 

= dx [\4> k (x)\ 2 - \<t)* k A^ A<f) k \ + surface term. (41) 
J u k 

The first term simply vanishes because of (40). The surface term results from the fact that 
in the above expression the operator A^ is only formally the hermitian conjugate of A. The 
difference of the spectral densities is therefore given by 

1 C 2m 

-9\k) = -^Jdx d x (V\Mz)\ 2 ) = -p^-f ( 42 ) 

where we have omitted a term <f)*d x <j) because it is is even in x for large \x\ and thus only 
needed the asymptotic values U'[ipk(x = ±00)] = U'(±v) = ±m and that the mode func- 
tions are plane waves asymptotically, i.e. \<pk\ 2 1- This result coincides with the one from 
a direct evaluation of (39). 

Thus we obtain for the remaining correction from (37) 

am (1) 1 r dkd d -H£\, /7 , 2 r(^) , 

2 



L d - X 4 J (27r) d u v ' d(4,r) 



This reproduces the correct, known result for the susy kink mass correction 
—m/(2rc) (for d — 1) and the surface (string) tension of the 2+1 dimensional susy kink 
domain wall AM^/L = -m 2 /(8n) (for d = 2) [21]. 

Notice that the entire result is produced by an integrand proportional to the extra mo- 
mentum component £ 2 , which for strictly d — 1 would not exist. 



10 The zero mode contributions in fact do not cancel by themselves between bosons and fermions, because 
the latter are chiral. This non-cancellation is in fact crucial in energy cutoff regularization (see Ref. [21]). 
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E. Anomalous contributions to the central charge and extra momentum 

In a kink (domain wall) background with only nontrivial x dependence, the central charge 
density Z x receives nontrivial contributions. Expanding Z x around the kink background 
gives 

Z x = Ud x <p K - ^=d x <p K + d x (U V ) + Id^U'rf) + 0{rf). (44) 
v 2A * 

Again only the part quadratic in the fluctuations contributes to the integrated quantity 
at one-loop order 11 . However, this leads just to the contribution shown in (38), which 
matches precisely the counterterm mSv 2 from requiring vanishing tadpoles. Straightforward 
application of the rules of dimensional regularization thus leads to a null result for the net 
one-loop correction to (Z x ) in the same way as found in Refs. [5, 6, 8] in other schemes. 

On the other hand, by considering the less singular combination (H+Z x ) and showing that 
it vanishes exactly, it was concluded in Ref. [9] that (Z x ) has to compensate any nontrivial 
result for (if), which in Ref. [9] was obtained by subtracting successive Born approximations 
for scattering phase shifts. In fact, Ref. [9] explicitly demonstrates how to rewrite (Z x ) into 
— (if), apparently without the need for the anomalous terms in the quantum central charge 
operator postulated in Ref. [11]. 

Because the authors of Ref. [9] did not discuss regularization of (Z x ), the manipulations 
needed to rewrite it as —(H) (which eventually is regularized and renormalized) are not 
defined in their work. If we choose to use dimensional regularization, (Z x ) contains the 
mode energies oo = \Jk 2 + m 2 + £ 2 instead of uot = Vk 2 + m 2 and so the manipulations 
carried through in Ref. [9] (eq. 56) are no longer possible. Using dimensional regularization 
one in fact obtains a nonzero result for (H + Z x ), apparently in violation of susy. However, 
dimensional regularization by embedding the kink as a domain wall in (up to) one higher 
dimension, which preserves susy, instead leads to 

(H + Z x -P y ) = 0, (45) 

i.e. the saturation of (24), as we shall now verify. 

The bosonic contribution to (P y ) involves an I integral which is scale-less and odd and 
thus vanishes. Only the fermions turn out to give interesting contributions: 

(Vy) = | «.»<)„</.} 

^/(S^^'g(>^l^) + ^W 2 -N 2 )). (46) 

We have already omitted the contributions which vanish either by symmetric integration 
or due to scale-less integrals, which are zero in dimensional regularization. The remaining 
^-integration no longer factorizes because oo = \/k 2 + I 2 + m 2 , and is in fact identical to the 
finite contribution in (Ti) obtained already in (37): 

f \ 1 f dkd d ~He 2 nU ,, 2T(^) , 

-az = J dx { v,) = - J -^r-m = - d -^ m -K (47) 



Again, this does not hold for the central charge density locally [11, 12]. 
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So for all d < 2 we have BPS saturation, (H) = | (Z x — P y )\, which in the limit d — > 1, the 
susy kink, is made possible by a non- vanishing (-P^). The anomaly in the central charge is 
seen to arise from a parity- violating contribution in d — 1 + e dimensions which is the price 
to be paid for preserving supersymmetry when going up in dimensions to embed the susy 
kink as a domain wall. 

It is again the difference in the spectral densities, 9', which determines the one-loop 
corrections, which thus depend only on the derivative of the pre-potential (or equivalently 
the second derivative of super potential W — J dipU (</?)) at the critical points ±v. In general 
the spectral density difference for a model with spontaneously broken Z 2 symmetry is given 
by 

W'(v) - W"(-v) 

'<*> - v^ww ' (48) 

which has an obvious generalization for Zjv symmetric models like the sine-Gordon model. 
From (40,41) one can see that this quantity is closely related to the index of the operator 
AA^ . For the simple models considered here, where only one spatial direction is nontrivial, 
9'(k) is easily obtained from the Dirac equation in the asymptotic regions x — > ±oo, far 
away from the kink [6]. But as we will see below, in case of a less trivial background like 
the vortex, the formulation as surface term will provide essential simplifications [31]. 



F. Dimensional reduction and evanescent counterterms 

In the above, we have effectively used the 't Hooft-Veltman version of dimensional regu- 
larization [32] in which the space-time dimensionality n is made larger than the dimension 
of interest. This is possible in a supersymmetric way if the model of interest can be obtained 
from a higher dimensional supersymmetric model by dimensional reduction. The nontrivial 
corrections to the central charge of the kink come from the "genuine" momentum operator 
P y , and are due to a spontaneous breaking of parity. 

We now comment on how the central charge anomaly can be recovered from Siegel's 
version of dimensional regularization [33-35] where n is smaller than the dimension of space- 
time and where one keeps the number of field components fixed, but lowers the number 
of coordinates and momenta from 2 to n < 2. At the one- loop level one encounters 2- 
dimensional b 1 ^ coming from Dirac matrices, and n-dimensional 5^ from loop momenta. 
An important concept which is going to play a role are the evanescent counterterms [36-38] 

involving the factor ^S^y^ , where 5^ = — has only e = 2 — n nonvanishing components. 

The supercurrent is given by j M = — (flip + [/(</?)) 7^. In the trivial vacuum, expanding 
into quantum fields yields 

j, = -(&, + U'(v) V + ^"{v) rf^j ^ + ^=^ 2 7^- ( 49 ) 

Only matrix elements with one external fermion are divergent. The term involving U"(v)r) 2 
in (49) gives rise to a divergent scalar tadpole that is cancelled completely by the counterterm 
<5/i 2 (which itself is due to an rj and a ip loop). The only other divergent diagram is due to 
the term involving fir] in (49) and has the form of a ip- self energy. Its singular part reads 

(0b» div = iU"(v) f dx [ f" 2 f?"\ M 212 "(P)- (50) 

J Q J (2-n) n \k 1 + p z x(l — x) + m 2 \ 2 
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Using S^ = 5^ — 5^ we find that under the integral 

m\ P f v = ^ti^(p). (si) 

Hence, the regularized one-loop contribution to the susy current contains the evanescent 
operator 



so that 



(52, 



It is called evanescent because the numerator vanishes in strictly n = 2; for n ^ 2 it has 
a pole, but in matrix elements this composite operator gives a finite contribution. j^ lv is 
by itself a conserved quantity, because all fields depend only on the n-dimensional coordi- 
nates, but it has a nonvanishing contraction with 7 M . The latter gives rise to an anomalous 
contribution to the renormalized conformal-susy current $j™ n ' where j^ en ' = — j^ lv , 

^(^TOanom. = -7 M J^ = ~1>. (53) 

(There are also nonvanishing nonanomalous contributions to d^(^j^) because our model is 
not conformal-susy invariant at the classical level [39, 40].) 

Ordinary susy on the other hand is unbroken; there is no anomaly in the divergence 
of J™ 1 '. A susy variation of j M involves the energy-momentum tensor and the topological 
central-charge current according to 

5j M = -2T/ 7 ^-2^7 5 e, (54) 

where classically = e^ u Ud u ip. 

At the quantum level, the counter-term = — j^ v - induces an additional contribution 
to the central charge current 

c iom = l~ e ^ dPU ' ( 55 ) 

which despite appearances is a finite quantity: using that total antisymmetrization of the 
three lower indices has to vanish in two dimensions gives 

and together with the fact the U' only depends on n-dimensional coordinates this finally 
yields 

C™ = (57) 

in agreement with the anomaly in the central charge as obtained previously. 

We emphasize that itself does not require the subtraction of an evanescent counterterm. 
The latter only appears in the susy current j^, which gives rise to a conformal-susy anomaly 
in /tj^. A susy variation of the latter shows that it forms a conformal current multiplet 
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involving besides the dilatation current T^x" and the Lorentz current T p u x p e vp also a current 
jfj) 1 = xPe pi^C fl - We identify this with the conformal central-charge current, which is to be 
distinguished from the ordinary central-charge current 

Since d^Pj 1 = e^C M 5 an d e M1 , is invertible, the entire central-charge current £ M enters in 
the divergence of the conformal central-charge current, whereas in the case of the conformal- 
susy current it was the contraction 7 M j M . 

The current thus has the curious property of being completely determined by its own 
divergence. For this reason it is in fact not associated with any continuous symmetry (as is 
also the case for the ordinary central-charge current, which is of topological origin). In the 
absence of classical breaking of conformal invariance it is conserved trivially by its complete 
disappearance and then there is no symmetry generating charge operator. Nevertheless, in 
the conformally noninvariant susy kink model this current is nonvanishing and has in addi- 
tion to its nonanomalous divergence an anomalous one, namely the anomalous contribution 
to the central charge current inherited from the evanescent counterterm in the renormalized 
susy current. 



G. Multiplet shortening, BPS saturation and fermion parity 

We construct representations of the strictly two dimensional algebra 

Q +2 = H + Z , Q' 2 = H-Z , {Q + ,Q-} = 2P, (58) 

where Q ±J< = Q d= are hermitian, by going to the rest frame (P = 0) in the topological 
sector, i.e. with non-vanishing central charge. M and Z are then ordinary numbers. In the 
general case, M ^ \Z\, the irreducible representations are two-dimensional {|E_)|,E+)}. 
The supercharges can be represented as 

Q + = VM + Z (n , Q~ = VM-Z a 2 . (59) 

For BPS states the absolute value of the central charge is per definition equal to the 
energy of this state, i.e. for the eigen values in (58) we have 

M-\Z\=0. (60) 

We choose Z = — M which corresponds in our convention to the kink. The algebra (58) in 
such a BPS representation becomes now 

Q +2 = , Q- 2 = 2M , {Q+,Q-} = 0. (61) 

Because of the hermiticity of Q + , (61) implies ||<5 + |S)|| 2 = and thus 

g + |S) = 0. (62) 

This equation is equivalent to (60) for the definition of BPS states and means that the BPS 
states are left invariant by half of the supersymmetry, namely Q + in our case. Operators 
and states can be characterized by the cohomology of the operator Q + . Analogous to BRST 
exact operators which have vanishing matrix elements for physical states we can say that 
each operator which is Q + - exact has vanishing expectation value for BPS states: 



O = {Q + , O'} (BPS\0\BPS) = 0. 



(63) 
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For the other supercharge in (61) Q , which acts nontrivially, there exist two inequivalent 
irreducible representations, 

Q-\E) = ±V2M\E), (64) 

which are connected by a Z 2 transformation ip — > — xjj for all fermions, which is clearly a sym- 
metry for each fermionic action. Thus the irreducible representation is one-dimensional and 
the fermionic operator is diagonal [41]. This is the reason why it was originally thought that 
multiplet shortening does not occur in two dimensions [4, 6, 11]. Therefore a reducible two 
dimensional representation for the soliton states was assumed such that the fermion parity 
operator (— 1) F is still defined. For a reducible two dimensional representation {|Xb), |£f)} 
we may choose: 

Q- = V2Ma x , (-l) F = a 3 , (65) 

so that (— 1) F is diagonal in this representation and Q~ has fermion parity —1, i.e. 
{Q~,(-1) F } = 0. Note that this is the direct sum of the two inequivalent irreducible 
representations (64), which are obtained as |(|S b ) ± |E f )). 

Witten and Olive [4] argued that in four dimensional susy gauge theories the number 
of particle states is not changed in the Higgs phase, although massive representations have 
2^ times as many states than massless one. Thus, they concluded that the Higgs phase 
corresponds to a BPS saturated representation which has the same number of physical states 
as the massless representation. Because of this multiplet shortening the BPS saturation 
should be protected against perturbative corrections since they cannot change the number 
of particle states. 

The counting of susy soliton states in two dimensions is somewhat peculiar (see below) 
and the loss of fermion parity (64) suggested a two dimensional representation, as for the 
non-susy soliton [18], and thus no multiplet shortening would occur. In [11], nevertheless 
BPS-saturation was assumed, to match the central charge correction to the mass correction 
obtained in [8]. The crucial relation for BPS saturation is the annihilation by one super 
charge (62). It was stated that this relation is protected without multiplet shortening, by 
analogous arguments that constrain supersymmetry breaking [29] . A simple argument shows 
that this is not sufficient. Assume that in some approximation a reducible multiplet is BPS 
saturated, i.e. Q + \T,i) = 0. Since the operator Q + is hermitian its representation is of the 
form (because of the reducibility the Z 2 -grading through (— 1) F exists) 

and the BPS states can be separated in zero-eigen states of M and M\ To answer the ques- 
tion if the multiplet remains BPS saturated under perturbations (corrections) we consider 
the quantity 12 



lim Tr 



-l) F e- 13 Q+2 =n° h -n° { = Ind(M), (67) 



which is the difference between the number of zero-eigenvalue eigen states (singlets) of M 
and M\ That this index is invariant under perturbative corrections can be seen quite 
analogously to the arguments of [29]. If a state is no longer annihilated by Q + , say for 



12 Note that BPS states (Q +2 = 0) contribute Tr\ B ps(-l) F = n° — n° f , whereas for non-BPS states is 
Q +2 = (H + Z) > such that their contribution vanish for (3 — > oo. 
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example one with fermionic parity, then also the bosonic super-partner state is no longer 
annihilated: 

+ Q+|f > ~ Q + Q~\b) = -Q-Q + \b) -> Q + \b) + 0. (68) 

So the difference in the number of singlet states is unchanged under perturbative corrections, 
and thus it can be calculated in a semi-classical approximation. So what can this index tell 
us? In case that it would be nonzero, there would exist, at least one, BPS saturated state, 
which is then protected against quantum corrections. But in this case (— 1) F is no longer 
defined as we will see immediately. If the index vanishes in an approximation, n° — n° = 0, 
the number of fermionic and bosonic singlets coincides. The trivial case is of course that they 
both vanish already in the approximation and there are no BPS states. In the nontrivial case 
there exist susy pairs of BPS singlets in the approximation, but susy does not protect them 
from being lifted pairwise above the BPS bound as described in (68). But this is exactly 
the case of the Z 2 symmetric two-dimensional multiplet (65). So the equality between the 
mass correction and the anomalous contribution to the central charge needs a different 
explanation. In fact, it was found that one has to give up the usual fermion parity for 
the topological soliton state which is then a single-state short super-multiplet (64) [41]. If 
we look now again on the BPS saturation equation (62), we see immediately that a lift 
above the BPS bound would give a twice as long irreducible multiplet (59) which cannot be 
caused by perturbative corrections. So, in the absence of other mechanisms as for example 
a difference in a conserved quantum number, multiplet shortening is a necessary condition 
for BPS saturation being protected. 

Up to now we have only discussed abstract representations of the susy algebra (58). In a 
quantum field theory the operators in the algebra (58) and their representations correspond 
to Heisenberg operators of conserved, i.e. time- independent, charges and Heisenberg states. 
In general neither operators nor states in the Heisenberg picture are known explicitly. Instead 
one quantizes the field operators in the interaction picture in terms of creation/annihilation 
operators which are defined w.r.t. a perturbative ground state. The canonical commutation 
relations imply an algebra for the mode coefficients which usually has to be represented in an 
irreducible manner. This determines which kind of the above representations is realized in 
the (perturbative) quantum field theory. From the fermionic creation/annihilation operators 
(34) one obtains an infinite dimensional Clifford algebra of pairs of generators 7^ = (bk + b k ) 
and 7^ = i(b k — b k ) and the single generator b corresponding to the zero-mode. This is 
quite analogous to an odd dimensional Clifford algebra, and the operator b plays the role 
of the 75 of the even-dimensional algebra 7^. So there are two inequivalent representations 
of the full algebra governed by the sign of the "gamma five" operator 6 - Because bo has to 
anti-commute with the fr^'s it cannot be represented as a number, as in the quasi-classical 
approximation [41]. The &fc-algebra can b e represented as usual on a Fock space, constructed 
from the Clifford vacuum \Q) with bk\ty = 0. The whole algebra, including b can then be 
realized by two inequivalent irreducible representations [14]: 

|s±) = i(l±6o)|fi) , b \s ± ) = ±\s ± ) , 6 fc |s±> = 0. (69) 

According to the usual fermion-parity counting 6 is an odd, i.e. fermionic operator, and 
thus the ground states \s±) are half fermionic and half bosonic. But in two dimensions there 
is less distinction between fermions and bosons. In fact, since there are no rotations in one 
spatial dimension, the definition of fermion parity is more abstract. Indeed, the assignment 
of fermion number to different vacua depends on the sign of the eigenvalue of the fermion 
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mass matrix at the considered vacuum [29]. In the case of the kink this means that if the 
vacuum with tp = +v is defined to be bosonic the vacuum at — v is automatically fermionic. 
Now a topological state like the kink connects these two vacua with opposite fermion parity 
which heuristically explains that this state cannot have a definite fermion parity in the usual 
sense. 

Let us now check if also semi-classically the BPS saturation condition is satisfied. With 
the regularized mode expansion for the quantum fields (30,32) and the BPS equation d x ipx + 
U = one obtains 

Q + \s±) = [dx [^-dy^ + {d x r 1 + U'r 1 )4,~]\s ± ) + 0{h) 

AU 

— (V^i-V^+i)albl\s ± ) = 0. (70) 
In 

So both states are BPS saturated semi-classically. 



_ . j d a - L £ ^ 
= 7 (2^4P 



IV. SUPERSYMMETRIC VORTICES IN 2+1 DIMENSIONS 

Following [42] we next consider the Abrikosov-Nielsen-Olesen [43-46] vortex solution of 
the abelian Higgs model in 2+1 dimensions which has a supersymmetric extension [47, 
48] (for related models see [49, 50]) such that classically the Bogomolnyi bound [51] is 
saturated. We implement dimensional regularization of the 2 + 1-dimensional M = 2 vortex 
by dimensionally reducing the M = 1 abelian Higgs model in 3 + 1 dimensions. We confirm 
the results of [47, 52, 53] that in a particular gauge (background-covariant Feynman-'t 
Hooft) the sums over zero-point energies of fluctuations in the vortex background cancel 
completely, but contrary to [47, 52] we find a nonvanishing quantum correction to the vortex 
mass coming from a finite renormalization of the expectation value of the Higgs field in this 
gauge [31, 53]. In contrast to [47], where a null result for the quantum corrections to the 
central charge was stated, we show that the central charge receives also a net nonvanishing 
quantum correction, namely from a nontrivial phase in the fluctuations of the Higgs field in 
the vortex background, which contributes to the central charge even though the latter is a 
surface term that can be evaluated far away from the vortex. The correction to the central 
charge exactly matches the correction to the mass of the vortex. 

In Ref. [52], it was claimed that the usual multiplet shortening arguments which prove 
BPS saturation would not be applicable to the Af = 2 vortex since in the vortex background 
there would be two rather than one fermionic zero modes [54], leading to two short multiplets 
which have the same number of states as one long multiplet. 13 We show however that the 
extra zero mode postulated in [52] has to be discarded because its gaugino component is 
singular, and that only after doing so there is agreement with the results from index theorems 
[54-56]. For this reason, standard multiplet shortening arguments do apply, explaining the 
BPS saturation at the quantum level that we observe in our explicit one-loop calculations. 



Incidentally, Refs. [52, 54] considered the supersymmetric abelian Higgs model augmented by a Chern- 
Simons term. 
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A. The model 

The superspace action for the vortex in terms of 3+1-dimensional superfields contains an 
Af = 1 abelian vector multiplet and an Af = 1 scalar multiplet, coupled as usual, together 
with a Fayet-Iliopoulos term but without superpotential, 

C = J d 2 9W a W a + J d 4 9$e- eV $ + k J d 4 9V. (71) 

In terms of 2-component spinors in 3+1 dimensions, the action reads 14 

£ = -\n» + ^-<^x P + \D* + {K-e\^)D 

-\D^\ 2 + riolpD^ + \F\ 2 + V2e [<P*x«r + <l>xd & ] , (72) 

where = <9 M — ieA^ when acting on <fi and if, and = d fl A u — dyA^. Elimination of 
the auxiliary field D yields the scalar potential V = \D 2 = |e 2 (|0| 2 — v 2 ) 2 with v 2 = n/e. 
This model is invariant under the following transformation rules: 

SAp = e a af X a ~ ^ aX a , SD = i(e a a^% X a + ^^Ax") (73) 
5 X a = -iF^ va ^ - De a , Sx& = iF IJH/ a' w Je $ -De & , (74) 

for the gauge multiplet; the matter multiplet transforms as 

50 = -V2e^ a , S^ a = -iyft&a^D^ + V2F*e a 
5<f* = V2^ & , 5$ & = -iV2(D^Ya» a t a + V2Ft^ 

5F = -iV2e a a^D^ & + 2e<f*e aX a , 6F* = -iy/2^ a D^ a + 2e0e^, (75) 

where 

a^p = \(^% ~ ° vaP o%) , = - (76) 

As always, before eliminating auxiliary fields, the transformation rules of the gauge multiplet 
do not depend on the matter fields, and the transformation rules themselves are gauge 
covariant and lead to a closed superalgebra. 

In 2+1 dimensions, dimensional reduction gives an Af = 2 model involving, in the notation 
of [52], a real scalar N = A 3 and two Dirac spinors if = (-0°), x — (x a )- 

Completing squares in the bosonic part of the classical Hamiltonian density for time- 
independent fields one finds the Bogomolnyi equations and the central charge 

H = \F 2 k i + \D k 4>? + \e 2 {\<t>\ 2 -v 2 Y 

1 \D k <f + ie kl D x <f>\ 2 + \ (F 12 + e(|0| 2 - v 2 )f 



2' 



+^v 2 e kl F kl - id k {e kl (t>*D l( t>) with k,l — l, 2. (77) 



14 



Our conventions are rf v = (—1, +1, +1, +1), x a = e Q/3 X/3 an d X° — e<x ^Xp with e al3 = e a p = -e Qf9 = — e & p 
and e 12 = +1. In particular we have ffa — (ipa)* but fp a — — (ip a )*. Furthermore, <7^ = (— 1,<?) with 
the usual representation for the Pauli matrices a, and a^" 13 — cr^ 130 with = (!,<?)• 
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The classical central charge reads 

Z = J d 2 xe kl d k (ev 2 A l - ipDrf) , (78) 

where asymptotically Di<p tends to zero exponentially fast. Classically, BPS saturation 
E — \Z\ — 2nv 2 \n\ holds when the BPS equations (Di ± iD 2 )(f) = D±(f> = and F 12 ± 
e(\(f)\ 2 — v 2 ) = are satisfied, where the upper and lower sign corresponds to vortex and 
antivortex, respectively. In this paper we use the vortex solution. For winding number n, it 
is given by 

(j)v = e me f(r), eA Y + = -ie te ^±^, J^ = A* ±iJ% (79) 

or, alternatively, 

*, = rt + i tf = p±^)" /(r) , e^ = !^^, (80) 

where f'(r) = ff(r) and a'{r) = re 2 (/(r) 2 — v 2 ) with boundary conditions [46] 

a(r — > oo) = 0, f(r — > oo) = v, 
a(r -> 0) = n + 0(r 2 ), /(r -> 0) -> r n + 0(r"+ 2 ). (81) 



B. Fluctuation equations 

For the calculation of quantum corrections to a vortex solution we decompose into a 
classical background part 0v and a quantum part 77. Similarly, A^ is decomposed as A^ + a^, 
where only with /x = 1, 2 is nonvanishing. We use a background i?£ gauge fixing term 
[57, 58] which is quadratic in the quantum fields, 

£g.fix = - ie£(<fori* - <PwV)) 2 - (82) 

The corresponding Faddeev-Popov Lagrangian reads 

£ g host = b {&l - e 2 ^ {2 |0 V | 2 + <PW + <PvV})c (83) 

The fluctuation equations in 2+1 dimensions have been given in [47, 52] for the choice 
£ = 1 (Feynman-'t Hooft gauge) which leads to important simplifications. We shall mostly 
use this gauge choice when considering fluctuations in the solitonic background, but will 
carry out renormalization in the trivial vacuum for general £ to highlight some of the gauge 
dependences. 

Because we are going to consider dimensional regularization by dimensional reduction 
of the 3+1 dimensional model, we shall need the form of the fluctuation equations with 
derivatives in the x 3 direction included. (This one trivial extra dimension will eventually be 
turned into e — > dimensions.) 

In the 't Hooft-Feynman gauge, the part of the bosonic action quadratic in the quantum 
fields reads (m = 0, 1, 2 but pi — 0, 1, 2, 3) 

411 = ~\{d,a m ) 2 - \{d,a,) 2 - e 2 \M 2 al 

-\D^\ 2 - e 2 (3|0 v | 2 - v 2 )\ V \ 2 - 2iea» [rj*D^ y - r,(D^th) m ] ■ (84) 
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In the trivial vacuum, which corresponds to 0v — *■ v an d ^ — > 0, the last term vanishes, 
but in the solitonic vacuum it couples the linearized field equations for the fluctuations 
B = (77, a + /\/2) with a+ = a\ + ia-i to each other according to (k = 1, 2) 

^- d t) B -{ -iyfaiD-W -d 2 + 2e 2 |0 v | 2 J { ) 

The quartet (03, ao, b, c) with b, c the Faddeev-Popov ghost fields has diagonal field equations 
at the linearized level 

(dj - 2e 2 |0v| 2 )Q = 0, Q = (a 3 ,a ,b,c). (86) 

For the fermionic fluctuations, which we group as U = p^, V = (%), the linearized field 
equations read 

LU = i(d t + ds)V, tfV = i(d t - d 3 )U, (87) 

with 

-^y (88 ) 

-V2e0y id- J \ v2e0y l 9+ J 

Iteration shows that U satisfies the same second order equations as the bosonic fluctua- 
tions B, 

L^LU = (<9 3 2 - <9 2 )£7, iJLB = (<9 3 2 - df)B (89) 
WV = (dl - d?)V, (90) 

with L^L given by (85), whereas V is governed by a diagonal equation with 

rrt_f-W + ^v| 2 + eV \ 

V -dl + 2e^ v \ 2 )- { } 

(In deriving these fluctuation equations we used the BPS equations throughout.) 



C. Renormalization 

At the classical level, the energy and central charge of vortices are multiples of 2nv 2 with 
v 2 = k/c Renormalization of v 2 , even when only by finite amounts, will therefore contribute 
directly to the quantum mass and central charge of the M = 2 vortex, a fact that has been 
neglected in the original literature [47, 52] on quantum corrections to the M = 2 vortex. 

In 2+1 dimensions, it is possible, just as in the case of the 1+1-dimensional supersym- 
metric kink, to adopt a "minimal" renormalization scheme where the wave function renor- 
malization constants are set to unity, and only v 2 is renormalized. The renormalization of 
v 2 is then fixed by the requirement of vanishing tadpoles in the trivial sector of the 2+1 
dimensional model. The calculation can be conveniently performed by using dimensional 
regularization of the 3+1 dimensional M = 1 model. By going down in the number of spa- 
tial dimensions, 3 — > 3 — (1 — e), and setting e — > eventually, we have a supersymmetry 
preserving regularization method in analogy to the embedding of the supersymmetric kink 
as a domain wall in up to 2 + 1 dimensions. 
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For the calculation of the tadpoles we decompose = v + r\ = v + (cr + ip) / y/2, where 
o is the Higgs field and p the would-be Goldstone boson. With the gauge fixing term 
(82) and £ = 1 all fields (including the Faddeev-Popov fields b and c) have the same mass 
| to | = y/2ev. This gauge choice also avoids mixed a^-p propagators, but there are mixed 
X~ip propagators, which can be diagonalized 15 by introducing new spinors s — {ip + ix) / V% 
and d — (ip — ix)/V% with mass terms ^(s a s a — d a d a ) + h.c. 

The part of the interaction Lagrangian which is relevant for a tadpoles to one-loop order 
is given by 



r'mt 
a— tadpoles 



em 



= e( X *r + X^ a ) o - — (a 2 + p 1 ) cr - eTO(a 2 + the - 5v 2 ) a, (92) 



where b and c are the Faddeev-Popov fields. 

The one-loop contributions to the a tadpole thus read 




a 0, 1,2,3 





b,C 



5v z 



-em) x 



(93) 



{-2trl 2 /(m] 



TO 



+ - J(^to) + [3/(to) + ^l(^m) 



-ZI(C*m)-6v 2 }, 



where 



/(to) 



/ 



d^k 



TO 



1+e 



r(-l-l) 



(27r) 3+e k 2 + to 2 (4 7 r) 1 + e / 2 r(-|) 



TO 
47T 



+ 0(6). 



(94) 



Note that because we use dimensional reduction the component 03 is kept also in the limit 
of 2+1 dimensions. In fact, the result equals that of ordinary dimensional regularization in 
2+1 dimensions, where a 3 appears as an additional scalar field. 

Requiring that the sum of tadpole diagrams (93) vanishes fixes 5v 2 , 



Sv 2 



I(£*m] 



1+e 



D=3 



8tt 



-TO. 



(95) 



Since in dimensional regularization there are no poles in odd dimensions at the one-loop 
level, the result for Sv 2 is finite, but it is nonvanishing. Because the classical mass of the 
vortex is My = 2ttv 2 = irm 2 /e 2 , the counterterm 5v 2 is the only one that is of importance 
to the one- loop corrections to My. 

While the latter statement holds generally, the result for 5v 2 does depend on the renor- 
malization scheme adopted. If we consider more general renormalization schemes, which also 
allow for wave function renormalization, the latter do not affect the mass of the vortex, be- 
cause that is evaluated at a stationary point. However, they may affect the definition of 5v 2 . 
For example, Z^^l introduces a further counter-term to (93) such that the left-hand-side 



15 



The kinetic terms for s are —sj^d^s where s = s'a^ and 7^ = {— ia^, ia^a™} with m = 1,2,3, and 
likewise for d. Note also that {s^) a = —s a according to footnote 14. 
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of (95) is replaced by 5v 2 — v 2 5Z ( f l . While the residue of the scalar propagator is finite and 
does not enforce a nontrivial SZ^ to dispose of a UV divergence, one can nevertheless choose 
a finite value of SZ^ precisely such that Sv 2 = 0. By allowing for a nontrivial wave function 
renormalization of the vector boson field, Z 1 ^ 2 = Z' 1 , one can accommodate further renor- 
malization conditions (for example an on-shell definition of the elementary masses, residues, 
and coupling constants), which generically involve a nontrivial 5v 2 . In the following we shall 
however stick to the "minimal" scheme with Z^ — 1, which has been used predominantly also 
in the case of the 1+1 dimensional kink (for a detailed discussion of other renormalization 
schemes in the latter context see [21]). 

Since in general Sv 2 is gauge-parameter dependent, the remaining contributions to vortex 
mass and central charge must be gauge dependent, too, so that the final result is gauge 
independent when expressed in terms of physical parameters. (Notice that v 2 is not such a 
physical parameter, but in any well-defined renormalization scheme it can be related to the 
physical mass and coupling constant.) 

The gauge breaking term in fact breaks susy, but because the final result should be gauge- 
choice independent, we should get the correct result for the mass and central charge from 
this x-space susy-breaking approach. One could also use a superspace approach and fix the 
U(l) gauge symmetry without breaking rigid supersymmetry, but in the presence of solitons 
the background superspace formalism leads to some problems, as we sketch in the appendix. 

To gain further insight into the occurrence of tadpoles and their contribution to the 
quantum mass of the M = 2 vortex, we briefly consider an extension of the vortex model 
in which the U{\) anomaly in 3+1 dimensions cancels. In the 3+1 dimensional model, 
the U(l) coupling of the vector multiplet to the chiral multiplet is chiral (in terms of 4- 
component Majorana spinors it contains the matrix 75), and in order to cancel the chiral 
U(l) anomaly, additional scalar multiplets would be needed such that the sum over charges 
vanishes, X^ e « = 0- As it turns out, this anomaly is also of concern in 2+1 dimensions. 
There it does not entail a local gauge anomaly when £\ et 7^ 0, but leads to a parity- violating 
Cherns-Simon term [59-61]. 

The simplest extension fulfilling this condition is massless super-QED, consisting of a 
vector multiplet coupled to two chiral multiplets with opposite electric charges and still 
without superpotential [62]. In such a theory, the Fayet-Iliopoulos term is not generated in 
any order of perturbation theory if it is not present classically [63]. As we now briefly show, 
this does not imply the absence of quantum corrections to v 2 = nje if v 2 is nonzero. 

In massless super-QED, the kinetic term, Yukawa coupling and the scalar coupling to 
D for the new terms are the same as for and 0, except for an opposite overall sign. 
Eliminating D, the potential becomes 



where is the new complex scalar. We locate the vortex solution again in 0, and use 
the same gauge fixing term as in (82). 16 We find then two new couplings which produce 



In fact, it has been shown in [64] that the only solutions with nonvanishing winding number that are 
time-independent are the original vortex solutions, with <j> = and the vortex located in 4>. An analogous 
result holds also for the kink [65] . 




(96) 
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cr-tadpoles 

CX^,es = ^w(? 2 +^) (97) 

However, these new fields are massless and therefore the additional tadpole diagrams vanish 
in dimensional regularization. 

On the other hand, expanding about = the a field does not have tadpole diagrams- 
there is only a vertex left, but this does not give rise to a tadpole diagram because ip 
does not mix with x an d (or s and d) when <fi vanishes. 

Hence, the finite renormalization of v 2 in this extension of the 2+1-dimensional susy 
Higgs models is identical to that of the simpler model we have introduced above, provided 
the minimal renormalization scheme with trivial wave function renormalization constants is 
employed. 17 

In the following we shall show that the nonzero result for Sv 2 leads to a nonvanishing 
mass correction for the Af = 2 vortex in 2+1 dimensions and is in fact required to match 
an equally nonzero correction to the central charge in order that the BPS bound remains 
saturated. 

D. Quantum corrections to mass and central charge 

The expressions for the central charge and stress tensor can be constructed from the clas- 
sical action without any gauge artefacts. However, when one evaluates one-loop corrections, 
one uses the gauge-fixing term to obtain propagators and well-defined fluctuation equations, 
and then one should use the quantum expression for H . In that case one should consider 
sums over zero-point energies including unphysical degrees of freedom and Faddeev-Popov 
ghosts. This can be done in a well-defined manner by using dimensional regularization by di- 
mensional reduction from the 3+1 dimensional model. Using this method, the central charge 
contains the standard 2+1 dimensional terms and, as a potential anomalous contribution, a 
remainder from the momentum operator in the extra spatial dimension. 

1. Mass 

At the one-loop level, the quantum mass of a solitonic state is given by 

M = M cl + - ^bos - ^ Wferm + 6M ( 98 ) 

where M c \ is the classical mass expressed in terms of renormalized parameters, SM represents 
the effects of the counter-terms to these renormalized parameters, and the sums are over 
zero-point energies in the soliton background (the zero-point energies in the trivial vacuum, 
which one should subtract in principle, cancel in a susy theory). 



It should be noted, however, that the situation is more complicated in 3+1 dimensions, since there one 
cannot do without wave function renormalization. 
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In the £ = 1 gauge the sum over zero-point energies is formally 

= J^wt/ - X^ y ' (") 

where the quartet (03,00,6, c) cancels separately. (Note that in (99) all frequencies appear 
twice because all fields are complex.) 

Using dimensional regularization for models with solitons as developed in [21], these 
sums can be made well defined by replacing all eigen frequencies u>k in 2+1 dimensions by 
Wk,e — (k'jfe + ^ 2 ) 1//2 where i are the extra momenta, and integrating over I: 

5Z^~S^ = /(2^ / J^y Uk/ j rf Mk| 2 + K| 2 -M 2 -H 2 ](x;k) (100) 

where we have written out only the contributions from the continuous part of the spectrum, 
using that L^L and LL^ , which govern U and V, respectively, are isospectral up to zero 
modes. In dimensional regularization, the zero modes of L^L and LL^ become massless 
modes (with energy but they can be dropped because the £- integration is scaleless, 

and thus vanishes in dimensional regularization. 

The existing literature [47, 52] also proves that the spectral densities of the continuous 
spectrum is equal for U and V, so there is a complete cancellation of the sums over zero- 
point energies in the £ = 1 gauge. In analogy to the calculation performed in eqs. (41) and 
(42), we can verify this result explicitly as follows. 

Using that LU k = u>kVk and LW k = oukU k , we can write the difference in the spectral 
densities appearing in (100) as 

Ap(k) = J d 2 x [UlU k - U k V k ] = J d 2 x [UlU k - Lu k 2 (LU)lLU k ] (101) 

and partially integrate to obtain a surface term of the form 

Ap(k) = iuJ k~ 2 J d 2 x{d-[u\(iD + ui + v^e^v^)] + d+[u* 2 (id-U 2 - V^ecjA^ui)}} 

= iu7 2 lim & d9 [u\{-d e + in)u\ + u* 2 d e u 2 } = 0, (102) 

r— +00 J 

which vanishes because cyclindrical waves decay like |ii 1)2 | ~ r~ 1//2 . 

Hence, u u ~ Yl — 0, and all that remains is a possible renormalization of v 2 . As we 
discussed, in the minimal renormalization scheme where all wave functions renormalizations 
are trivial there is a finite result for 5v 2 , and this leads to a nonvanishing quantum correction 
of the vortex mass according to 



E = 2n\n\(v 2 + 5v% =1 ) = 2n\n\(v 2 - £) = \n\ - ^ 



(103) 



where £ = 1, since in other gauges the fluctuation equations for the B fields, i.e. r],a + , no 
longer match those of the U fermions. As always in quantum field theory, the explicit form of 
the quantum corrections depends on the definition of the parameters in the Lagrangian, i.e. 
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on the renormalization conditions employed, which we have chosen in the simplest possible 
manner. There are of course renormalization schemes which are more physical such as on- 
shell renormalization of the parameters in the trivial vacuum. In the case of kink and kink 
domain walls, an extensive analysis of renormalization schemes other than the minimal one 
can be found in [21]. 

The above result for the vortex mass in the minimal renormalization scheme agrees with 
[53] , where however a careful analysis of boundary conditions in the heat-kernel approach was 
needed because the vortex had to be put in a box to discretize the spectrum. In dimensional 
regularization one does not need to put the system in a box, and as a consequence there is 
no need to study the contributions from these artificial boundaries. 

In the super-QED model considered at the end of the previous section we have found that 
the additional (tilde) fields do not change 5v 2 as given in (95). Concerning the additional 
fluctuation equations for these fields, these are even simpler than that of the minimal model: 
The rj field is governed by 

(dl - d 2 )v = [-(Df) 2 + eV - |0v| 2 M= -D?D?V (104) 

where D* differs from D only in the sign in front of e. The fermionic tilde field equations 
read 

i£>*V = + W 2 , iD?il> 2 = i{d t -d 3 )il)\ (105) 
Iteration shows that rt> 2 has the same field equation as rj, so that in the mode sum we have 

J2u(D*jD*V) -J2u(D*VD*J). (106) 

The same arguments that prove Yl^u ~ Yl^v = ^2<^(L^L) — ^2u(LL^) = can now be 
repeated for these simpler operators, and so the above result for the mass correction to the 
M = 2 vortex in 2+1 dimensions remains unchanged when considering the vortex of the 
anomaly-free super-QED model. 



2. Central charge 

By starting from the susy algebra in 3+1 dimensions one can derive the central charge 
in 2+1 dimensions as the component T 03 of 

T ^ = - t -Tra^{Q t ,,r a } (107) 

where J v a is the susy Noether current. 

The antisymmetric part of T^ u gives the standard expression for the central charge den- 
sity, while the symmetric part is a genuine momentum in the extra dimension: 

(Z) = J d 2 x (T 03 ) = (z + P 3 ) • (108) 

(A similar decomposition is valid for the kink [22].) 

Z corresponds to the classical expression for the central charge. Being a surface term, its 
quantum corrections can be evaluated at infinity: 

(Z) = j d 2 xd k e kl {Q = jfW*)Uoo (109) 
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with (i = ev\Ai — i<fiDi<j) and v q = v 2 + 5v 2 . 

Expanding in quantum fields (j) — (f)y + rj, A — A y + a and using that the classical fields 
0v — ¥ ve md , Aj — > n/e, Djcpy — > as r — > oo, we obtain to one-loop order 



2tt 



(Z) = 2nnv 2 - 1 d9 <(<^ + rf)(Dj - tea e ){<f) Y + 77)) I™ 
./o 

/•27T 

= 27rn{t;o - (77*77) | r ->oo} - * / dO {(rfd e rj) - iecfAj (a e r]} - ie<p v (<W)} 

Jo 

= z a + z b (110) 

where we have used (r)(r— >oo)) — > (which determines <5t> 2 at infinity, but at finite r, the 
vacuum expectation value does not vanish [11, 12]), (a$) = 0, and {rfrjag) = 0(h 2 ). 

The first contribution, Z a , can be easily evaluated for arbitrary gauge parameter £, yield- 
ing 

Z a = 2itn{vl - i((a(x) + (p/o»| r — oo} 
= 27rnH 2 -i[/(m) + /(^m)]} 

= 2t™(^ - 5w 2 ) = 2ixnv 2 . (Ill) 

If this was all, the quantum corrections to Z would cancel, just as in the naive calculation 
of Z in the susy kink [5, 6]. 

The second contribution in (110), however, does not vanish when taking the limit r — > 00. 
This contribution is simplest in the £ = 1 gauge, where the 77 and ag fluctuations are governed 
by the fluctuation equations (85). In the limit r — > 00 one has |0v| —> v and -D_0v — > 
exponentially. This eliminates the contributions from (a e r]) (note that at finite r there is a 
cross-term in the kinetic terms for ag and 77). However, D 2 ., which governs the 77 fluctuations, 
contains long-range contributions from the vector potential. Making a separation of variables 
in r and 9 one finds that asymptotically 



\DU 2 ^\drV\ 2 + ^\(de-in)r)\* (112) 



so that the 77 fluctuations have an extra phase factor e md compared to those in the trivial 
vacuum. We thus have, in the £ = 1 gauge, 



—7 



2tt 



d9 {{rfdgrj) - ie(jf w (a e r]} - ie<f)y (<W)} 

= -7/ dO^de-n) £=1 = 27m (77*77) ^ = 2™^ . (113) 
Jo £ =1 

This is exactly the result for the one-loop correction to the mass of the vortex in eq. (103), 
implying saturation of the BPS bound provided that there are now no anomalous contri- 
butions to the central charge operator as there are in the case in the M = 1 susy kink 

[22]. ' 

In dimensional regularization by dimensional reduction from a higher-dimensional model 
such anomalous contributions to the central charge operator come from a finite remainder 
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of the extra momentum operator [22] . The latter is given by [52] 

Z c = (P 3 ) = J d 2 x(F 0t F 3t + (D <j>)iD 3 <j>+ (D 3 ^D <f> 

-iX°od 3 X ~ ii>o- D 3 ip) . (114) 

Inserting mode expansions for the quantum fields one immediately finds that the bosonic 
contributions vanish because of symmetry in the extra trivial dimension. However, this is 
not the case for the fermionic fields, which have a mode expansion of the form 



- 1 mttM^'* +d " * (115) 

where we have not written out explicitly the zero-modes (for which oo 2 = £ 2 ). The fermionic 



contribution to Z c reads 



Z c 



= /(^7^/^[KI 2 +KI 2 -hl 2 -|^| 2 ](^k) (116) 



where ui = + £ 2 , so that the i integral is nontrivial in dimensional regularization. Only 
the continuous spectrum can contribute because zero modes give scaleless integrals which 
vanish in dimensional regularization and if there were other discrete states, they would 
cancel between U and V. However, the x-integration over the mode functions and 
produces their spectral densities, and we find 

because Ap(k) = as we have seen in (102). Hence, \Z\ = \Z a + Z b \ = E, so that the BPS 
bound is saturated at the (one-loop) quantum level. 



E. Fermionic zero modes and multiplet shortening 

Massive representations of the Poincare supersymmetry algebra for which the absolute 
value of the central charge equals the energy, i.e. when the BPS bound is saturated, contain 
as many states as massless representations, which is half of that of massive representations 
for which the BPS bound is not saturated. These results also apply in 2+1 dimensions for 
the M = 2 super- Poincare algebra [52]. 

A particular multiplet of states is obtained by taking the vortex solution, and acting 
on it with the susy generators of the M = 2 susy algebra, which contains two complex 
charges Q + , Q~, and their hermitian conjugates (Q + Y and (Q~Y- One of these charges, 
Q + , annihilates the vortex solution, while the other one, Q~ , is to linear order in quantum 
operators proportional to the annihilation operator of a fermionic zero mode. 
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However, if there indeed is a second fermionic zero mode in the model as claimed in [52] 18 , 
in second quantization it would be present in the mode expansion of the fermionic quartet 
UandV, 



ai Xl + ai Xn + non-zero modes. (H8) 



As a result, there would then be a quartet of BPS states 

\v), a\\v), a\\v), a\a\\v) (119) 

comprising two short multiplets of M = 2 susy, which are degenerate and together have 
as many states as one long multiplet without BPS saturation. As stressed in [52], the 
standard argument for stability of BPS saturation under quantum corrections from multiplet 
shortening [4] thus would not be applicable. 

However, we shall now show that there is in fact only a single fermionic zero mode in a 
vortex background with winding number n = 1 [42]. To this end, we first observe that the 
zero modes must lie in U, because V is governed by the operator LV of Eq. (91), whose 
only zero mode solution is V = 0. A zero mode for U must satisfy LU = 0, and to analyse 
this equation we follow [52] and set ^(x,!/) = —ie l ^^^ +n ^ e u(r) and x 1 = e' l< - j+ ^ 9 d(r). The 
equation LU = reduces then to 

(d r -(a + j-l)/r V2ef \M 

V2ef d r + {j + \)/r)\d)-^ (120) 

where / = f(r) and a = a(r) satisfy /' = -/ and a' = re 2 (f 2 — v 2 ). Iterating this equation 
yields 

> + U r - ^£ - 2 e 2 / 2 ) 0. (121) 

Given a solution for u, the corresponding solution for d follows from LU = 0. 

For given j, this equation has two independent solutions, a linear combination of which 
yields solutions which decrease exponentially fast as r — > oo. Hence, both solutions should 
be regular at r = 0. For j ^ \, one has, using f(r — > 0) ~ r n , 

^ ~ u ~ r n {C x r^ + C 2 r-^-^) for r -> (122) 

which selects for n — 1 only j = —\. This solution is the zero mode that is obtained by 
acting with Q~ on the background solutions, which gives ip 1 = — «£L0 v /\/2 = —iy/2f, 
X 1 = F12 = — e(f 2 — v 2 ). For j = |, one finds for n — 1 near r = 

i) x ~ d (x + iy) + C 2 (x + iy)]nr . (123) 



In the literature one can in fact find two different conventions for indicating the number of fermionic zero 
modes. Like Refs. [52, 56] we only count the number of zero modes in the fermionic quartet (U, V) and 
not additionally those in the corresponding conjugated fields (U\V^). One zero mode in this way of 
counting then corresponds to a pair of creation/annihilation operators. Alternatively one may count the 
zero modes in both (U, V) and (U' , V') and thus ascribe one zero mode to each creation or annihilation 
operator. The latter way of counting is employed for instance in Ref. [50] . 
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For large r, ipi ~ e~ mr e te , as follows from (121). This solution corresponds to the second 
fermionic zero mode postulated in Ref. [52]. 

However, while (123) is finite at the origin, the associated gaugino component is not: 
(120) implies that 

f ~ C 2 — , (124) 

r 

so this solution has to be discarded when C 2 ^ 0. 

Similarly, one can show that for winding number n > 1 regularity of the gaugino com- 
ponent generically requires that j < — \ so that the correct quantization condition for 
normalizable fermionic zero modes is — n + \ < j < — |. Hence, there are n independent 
fermionic zero modes, not 2n as concluded in [52]. It is in fact only the former value that 
agrees with the results [54, 56] obtained from the index theorem [55]. 

As has been proved rigorously in [66], in the bosonic sector there are 2n zero modes, which 
are related to the above n independent fermionic zero modes by supersymmetry. In the R^=± 
background gauge d M a M — ie((f>yr]*— (j)yV) = 0, the bosonic zero modes satisfy a set of equations 
completely equivalent to those for the fermionic zero modes [54] . But the linearly dependent 
solutions (^) and correspond to linearly independent solutions for the bosonic zero 
modes a and rj. 19 In particular, for n — 1, the j = —\ solution (if) 1 = — ^(r),^; 1 = d(r)) 
with real u(r) and d(r) corresponds to the bosonic zero mode r](r) = —iu(r), (01,02) = 
(V2d(r), 0), while multiplying the fermionic solution by % corresponds to the bosonic zero 
mode r\[r) = u(r), (ai,a 2 ) = (0, v / 2<i(r)), which is evidently linearly independent of the 
former. For both solutions the R^=\ gauge condition is satisfied due to the lower component 
of the field equation (120). Conversely, one can start from the classical vortex solution and 
find two independent bosonic zero modes by considering their derivatives with respect to the 
x and y coordinates. Performing a gauge transformation to satisfy the R^=\ gauge condition 
leads one back to the above solutions. This additionally confirms that the above analysis 
has identified all fermionic zero modes in the quartet (U, V). 

We thus conclude that for the basic vortex (winding number n — 1) there is exactly one 
fermionic zero mode (corresponding to one pair of fermionic creation/annihilation operators) 
and this gives rise to a single short multiplet at the quantum level. Standard multiplet 
shortening arguments therefore do apply and explain the preservation of BPS saturation 
that we verified at one-loop order. 
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APPENDIX A: R ( GAUGE IN SUPERSPACE 

It seems natural to use superfields in the calculation of quantum corrections, since in 
this way we manifestly preserve rigid susy. Since we want to descirbe quantum fluctuations 
about a nontrivial background, we use the background formalism for superspace. The action 
reads i>e v ' ( I>, where $ = 4>\(x) +tp(x, 9) and V = Vy(x) +w(x, 9) for abelian gauge theories. 



For an analogous case see eq. (3.8) of Ref. [67]. 
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Expanding in terms of quantum fields, the classical action contains again off-diagonal terms 
quadratic in quantum fields, 



£kin. = J d 4 9 [(pyW^J + 1pW(f)y + ...]. 



(Af) 



To cancel these, we try to modify the usual D = 4 gauge fixing term D 2 VD 2 V into a 
superspace gauge fixing term 



£ fix = — (D 2 w + D 2 ^ct> Y ){D 2 w + D 2 ^ y ) 



(A2) 



In a trivial (constant) vacuum with eV = v, the terms ipcpy and (fiyip are antichiral and 
chiral, respectively. So to maintain these chirality properties we extend (f)y(x) in the soliton 
sector to an anti-chiial superfield (j>y(x, 6) = <py(x^ — iOa^O). Then 0y is chiral, and one has 
the identities 

i i 

(A3) 



D 2 D 2 ^<f> Y ) = ^ D 2 D 2 ^<j>y) = jfa. 



The terms — J d 4 9w(^(f)y + 4>(f)y) in the gauge-fixing term cancel the classical off-diagonal 
kinetic terms just as in x-space. 

However, one is left with the following diagonal terms in the gauge-fixing term 



/ 



d A e 



d i e$<fo-y><fr) = D 2 



{D 2 {U v ))^{Wy) 



(A4) 

Although D 2 D 2 (^(J)y) produces a factor □ which cancels the nonlocal !/□, further terms 
with one or two D a acting on ip(fi v yield nonlocal interactions. Thus it seems that one cannot 
construct a local gauge with £ = 1 in superspace. 20 
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